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Abstract 


SiiEplified  chemical-kinetic  mechanisms  were  employed  to  calculate  the  structure  of  the  dark 
zones  of  burning  double-base  and  nitramine  propellants.  These  reduced  mechanisms  are  expected 
to  be  useful  in  interior  ballistic  calculations.  First,  skeletal  mechanisms  comprised  of  22 
elementary  reactions  among  15  species  and  23  elementary  reactions  among  17  species  were  used 
to  calculate  dark  zone  structures  for  double-base  and  nitramine  propellants,  respectively.  The 
skeletal  mechanisms  were  previously  extracted  from  a  detailed  mechanism  incorporating  190 
elementary  reactions  involving  41  species.  Ignition  delay  times  (Xjg)  were  calculated  for 
homogeneous  mixtures  in  which  the  initial  concentrations  of  reactants  were  similar  to  those  found 
at  the  beginning  of  propellent  dark  zones.  The  tjg  calculations  were  performed  for  various  initial 
pressures  and  temperatures.  The  calculated  using  the  skeletal  and  detailed  mechanisms  agreed 
well.  Reduced  mechanisms  were  derived  from  the  skeletal  mechanisms  by  introducing  steady- 
state  approximations  for  a  number  of  species.  For  double-base  propellants,  a  reduced  mechanism 
of  three  global  reactions  was  obtained.  Reduced  mechanisms  utilizing  six  and  four  global 
reactions  were  deduced  for  nitramine  propellants.  The  and  structures  of  the  dark  zones 
calculated  using  the  reduced  three-step  mechanism  for  double-base  propellants  and  six-step 
mechanism  for  nitramine  propellants  were  in  agreement  with  calculations  made  using  the  skeletal 
mechanisms.  Agreement  using  the  four-step  reduced  mechanism  for  nitramine  propellants  was 
not  good;  however,  it  might  still  be  useful  for  some  apphcations. 


u 


ACKNOWLEDGMENTS 


The  research  at  the  University  of  California  at  San  Diego  was  supported  by  the  U.S.  Army 
Research  Office  through  grant  #ARO  DAAH-95-1-0108.  Dr.  W.  R.  Anderson’s  travel  to  present 
this  work  at  the  Twenty-Sixth  International  Symposium  on  Combustion  (July-August  1996)  was 
partially  supported  by  ERO  (Dr.  Roy  Reichenbach). 

The  authors  acknowledge  S.  W.  Bunte,  A.  M.  Dean,  J.  W.  Bozzelli,  R.  A.  Yetter,  F.  L.  Dryer, 
and  M.  C.  Lin  for  helpful  discussions  about  their  recent  research.  The  authors  are  also  grateful  to 
J.  W.  Bozzelli  for  providing  kinetic  estimates  for  the  HCN=HNC  reaction. 


Ill 


Intentionally  left  blank. 


IV 


TABLE  OF  CONTENTS 


Page 

ACKNOWLEDGMENTS  .  iii 

USTOFHGURES .  vii 

UST  OF  TABLES .  vii 

1.  INTRODUCTION . 1 

2.  SKELETAL  MECHANISMS  .  2 

3.  REDUCED  CHEMICAL-KINETIC  MECHANISMS .  6 

3.1  Double-Base  Propellents  .  8 

3.2  Nitramine  Propellents .  10 

4.  RESULTS  AND  DISCUSSION  .  12 

5.  SUMMARY  AND  CONCLUSIONS . 17 

6.  REFERENCES  .  19 

DISTRIBUTION  LIST .  21 

REPORT  DOCUMENTATION  PAGE .  25 


V 


fi^NTIONALLY  LEFT  BLANK. 


VI 


USTOFHGURES 


Figure  Page 

1 .  The  ratio  of  the  difference  between  the  rates  of  production  and  the  rates  of 
consumption  to  the  rates  of  production  for  the  species  N,  NH,  HNO,  HONO,  H, 

OH,  and  O  plotted  as  a  function  of  time .  9 

2.  The  ratio  of  the  difference  between  the  rates  of  production  and  the  rates  of 
consumption  to  the  rates  of  production  for  the  species  N2O  and  NO2  plotted  as 

a  function  of  time .  10 

3.  Profiles  of  temperature  (T)  and  mole  fractions  of  NO,  N2,  CO,  and  CO2  plotted 

as  a  function  of  time .  13 

4.  Profiles  of  temperature  (T)  and  mole  fractions  of  NO  and  N2  plotted  as  a  function 

of  time  .  14 

5.  Profile  of  temperature  (T)  and  mole  fractions  of  H2O  and  HCN  plotted  as  a 

function  of  time .  14 

UST  OF  TABLES 

Table  Page 

1 .  Skeletal  Chemical-Kinetic  Mechanisms  for  Modeling  the  Dark  Zone  of 

Double-Base  and  Nitramine  Propellents  and  Associated  Rate  Constants .  3 

2.  Comparison  of  the  Values  of  the  Ignition  Delay  Times  (tig)  Calculated  Using  the 

Various  Mechanisms  for  Double-Base  Propellents .  7 

3.  Comparison  of  the  Values  of  the  Ignition  Delay  Times  (Tjg)  Calculated  Using  the 

Various  Mechanisms  for  Nitramine  Propellents .  15 

4.  Comparison  of  the  Values  of  the  Ignition  Delay  Times  (Tjg)  Calculated  Using  the 

Various  Mechanisms  for  Nitramine  Propellents .  16 

vii 


Intentionally  left  blank. 


Vlll 


1.  INTRODUCTION 


Many  double-base  and  nitramine  propellents  exhibit  a  two-stage  flame  zone  during 
combustion.  In  this  process,  a  nonluminous  region  separates  the  primary  reaction  zone  near  the 
surface  of  the  propellant  from  the  luminous  secondary  flame  zone.  The  nonluminous  region  is 
commonly  referred  to  as  the  dark  zone.  In  the  primary  reaction  zone,  decomposition  of  the  solid 
propellent  takes  place,  and  moderately  reactive  intermediate  species  form.  These  intermediate 
species  convect  away  .from  the  surface  and  ignite  after  a  short  delay.  Recent  experiments  show 
that  chemical  effects  caused  by  the  use  of  propellents  of  differing  chemical  composition  can  lead 
to  major  differences  in  ignition  delay  times  (xQ  during  typical  large-caliber  gun  ballistic  cycles 
[1-4].  (The  Tjg  in  the  guns  are  not  the  same  as  the  dark-zone  chemical  delays.)  For  some 
propellents,  these  delays  are  undesirably  long.  It  is  believed  that  the  range  of  Xjg  for  various 
propellents  arises  from  differences  in  the  chemistry  taking  place  in  the  dark  zone  [3, 4]. 

It  has  been  established  that  the  dark  zone  of  solid  propellents  contains  large  amounts  of  NO, 
which  is  a  rather  weakly  reactive  oxidizer.  The  low  reactivity  of  NO  is  the  principal  reason  for 
the  dark-zone  formation.  This  weakly  reactive  mixture  is  rapidly  converted  to  equilibrium 
products  at  the  end  of  the  dark  zone.  Li  particular,  NO  is  converted  to  Nj,  and  the  temperature  at 
the  end  of  the  dark  zone  is  much  higher  than  at  the  beginning. 

The  temperature  and  species’  concentrations  change  very  slowly  with  distance  in  the  dark 
zone.  Sharp  gradients  in  these  profiles  are  observed  only  at  its  boundaries.  It  is  therefore 
reasonable  to  assume  that  the  structure  of  the  dark  zone  is  one-dimensional,  adiabatic,  and 
isobaric,  and  that  diffusion  of  heat  and  mass  are  negligibly  small. 

The  Tjg  for  any  reactive  mixture  is  defined  in  this  work  as  the  delay  time  from  the  initiation 
of  the  reaction  (time  zero)  to  the  time  of  maximum  heat  release  associated  with  the  thermal 
runaway.  It  is  possible,  in  unusual  situations,  to  obtain  an  incorrect  value  for  if  there  is  rapid 
heat  release  from  the  reaction  mixture  that  is  not  connected  with  the  steep  temperature  gradient 
generated  at  ignition.  All  x^g  reported  in  this  work  are  compared  with  temperature  profiles  to 
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ensure  validity  of  the  result.  Measured  dark-zone  thickness  is  directly  related  to  the  ignition 
delay  through  mass  conservation  and  gas  velocity  in  the  dark  zone  [5, 6]. 

A  detailed  chemical-kinetic  mechanism  has  been  developed  describing  the  time-dependent 
chemistry  for  reactive  mixtures  with  initial  compositions  at  temperature  and  pressure  conditions 
similar  to  those  at  the  beginning  of  the  dark  zone  of  double-base  and  nitramine  propellents  [7,  8] 
(the  first-stage  propellant  combustion  is  not  modeled).  These  reactive  systems  were  presumed  to 
be  homogeneous,  adiabatic,  and  isobaric.  Results  of  calculations  performed  using  this  detailed 
mechanism  have  also  been  compared  with  experimental  measurements  [7, 9]. 

In  this  work,  skeletal  mechanisms  comprising  fewer  elementary  reactions  than  those  in  the 
detailed  mechanisms  are  developed  for  use  in  modeling  the  structure  of  the  dark  zone  of 
propellents.  Reduced  chemical-kinetic  mechanisms  for  both  propellant  types  are  also  presented. 
Tjg  were  calculated  under  isobaric  and  adiabatic  conditions  using  the  skeletal  and  the  reduced 
mechanisms.  The  results  are  compared  with  those  calculated  using  the  detailed  mechanism. 
These  calculations  are  done  for  differing  values  of  pressure  (p)  and  initial  temperature  (T®). 
Pressure  was  held  constant  for  all  calculations. 

2.  SKELETAL  MECHANISMS 

The  skeletal  chemical-kinetic  mechanisms  employed  to  calculate  the  structure  of  the  dark 
zone  over  the  burning  surface  of  double-base  and  nitramine  propellents  are  shown  in  Table  1 . 

All  elementary  reactions  are  reversible.  In  Table  1,  the  elementary  reactions  1-22  among  15 
species  were  used  to  model  the  structure  of  the  dark  zone  of  double-base  propellents,  and 
reactions  1-12  and  23-33  among  17  species  were  used  to  model  the  stracture  of  the  dark  zone  of 
nitramine  propellents.  These  skeletal  mechanisms  were  derived  from  a  detailed  mechanism 
comprising  190  elementary  reactions  among  41  species  [9]  by  removing  those  reactions  found  to 
have  negligible  influence  on  the  calculated  values  of  the  Xjg.  The  skeletal  mechanisms  were  used 
in  the  development  of  the  reduced  mechanisms. 
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Table  1.  Skeletal  Chemical-Kinetic  Mechanisms  for  Modeling  the  Dark  Zone  of  Double-Base 
and  Nitramine  Propellents  and  Associated  Rate  Constants 


Note:  Units  are  kilocalories,  moles,  seconds,  and  cubic  centimeters.  Reactions  1-22  are  for  double-base 
propellents  ,and  reactions  1-12  and  23-33  are  for  nitramine  propellents. 
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Table  1.  Skeletal  Chemical-Kinetic  Mechanisms  for  Modeling  the  Dark  Zone  of  Double-Base 
and  Nitramine  Propellents  and  Associated  Rate  Constants  (Continued) 


No. 

Reaction 

An 

b„ 

En 

15. 

H2-i-N02=H0N0+H 

3.21E+12 

0.0 

28.81 

16. 

HN0+N02=H0N0+N0 

6.00E+11 

0.0 

1.987 

17. 

20H=0+H20 

6.00E+08 

1.3 

0.0 

18. 

H+0H+M=H20+M 
third-body  efficiency  /H2O  5.0  / 

1 

1.60E+22 

-2.0 

0.0 

19. 

N02+H=N0+0H 

1.30E+14 

0.0 

0.361 

20. 

N2O+ONO+NO 

6.92E+13 

0.0 

26.60 

21. 

N-hN0=N2+0 

3.27E+12 

0.3 

0.0 

22. 

NO+H=N-i-OH 

1.70E-1-14 

0.0 

48.80 

23. 

NH2+N0=N20-hH2 

5.00E-I-13 

0.0 

24.64 

24. 

C0+0+M=C02+M 

third-body  coefficiencies  /CO  1.77  /CO2  2.7  /H2O  5.0  /N2O 
5.0/ 

1 

2.36E-I-15 

0.0 

4.34 

25. 

HCN-l-0=NH+CO 

3.45E+03 

2.64 

4.98 

26. 

NH2+H=NH+H2 

4.00E+13 

0.0 

3.65 

27. 

NH2+NO-N2+H+OH 

9.30E+11 

0.0 

0.0 

28. 

NH2+N0=N2+H20 

2.00E-I-20 

-2.6 

0.924 

29. 

HNC+0=NH-hCO 

5.44E+12 

0.0 

0.0 

30. 

HNC+OH=HNCO+H 

2.80E+13 

0.0 

3.696 

31. 

HNC0+H=NH2+C0 

2.25E+07 

1.7 

3.80 

32. 

C0+N02=N0-hC02 

9.04E+13 

0.0 

33.78 

33. 

HCN-hM=HNC+M 

4.36E+26 

-3.34 

50.194 

Note:  Units  are  kilocalories,  moles,  seconds,  and  cubic  centimeters.  Reactions  1-22  are  for  double-base 
propellents  ,and  reactions  1-12  and  23-33  are  for  nitramine  propellents. 
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For  each  elementary  reaction  (n).  Table  1  presents  the  values  of  the  ■frequency  factor  (A„), 
the  temperature  exponent  (b„),  and  the  activation  energy  (EJ  that  appear  in  the  parametric 

b  A  -V 

expression  "exp[-Ejj/(RT)]  for  the  specific  reaction  rate  constant  (k„)  where  R 

represents  the  universal  gas  constant  and  T,  the  gas  temperature.  The  symbol  M  (in  reactions  1, 
12, 13, 14, 18, 24,  and  33)  represents  any  third  body,  and  the  catalytic  efficiencies  (qi)  of  various 
species  (i)  acting  as  the  third  body  in  these  reactions  are  shown.  The  catalytic  efficiency  of  those 
species  omitted  in  Table  1  is  assumed  to  be  unity.  The  concentration  of  the  third  body,  (Cj^)  is 

calculated  from  the  expression  Cj^  =  [p/(RT)]Sjlj  XjT^j  ,  where  Xj  is  the  mole  fraction  of  species 
i.  The  rates  of  the  elementary  reactions  1, 12, 13,  and  14  depend  on  pressure  and  are  calculated 
using  the  expression  [10]  k  =  F*k„-kL,  where  kL  -  +  koC^),  logioF  >=  logioFc/  <1  + 

[logio(koCM/k.„)]^},  Fc  =  ao  +  ajT.  The  values  of  k»,  ko,  a^,  and  aj  are  shown  in  Table  1. 

The  dark  zones  of  typical  double-base  propellents  are  known  to  include  large  amounts  of  NO, 
CO,  Hj,  Nj,  HjO,  and  COj  [6, 1 1].  Traces  of  CH4  and  C2H4  have  also  been  observed.  In  addition 
to  these  species,  large  amounts  of  HCN  and  trace  amounts  of  NjO  are  also  observed  in  the  dark 
zones  of  common  nitramine  propellents  [6].  Calculations  with  the  detailed  chemical-kinetic 
mechanisms  have  shown  that  the  trace  amounts  of  CH4  and  C2H4  have  negligible  influence  on  the 
Tjg  [8].  Inclusion  of  their  chemistry  would  unnecessarily  complicate  the  skeletal  mechanisms. 
Therefore,  these  species  are  omitted  from  the  initial  mixture.  The  main  difference  between  die 
skeletal  mechanism  for  double-base  and  nitramine  propellents  is  the  inclusion  of  HCN  chemistry 
in  the  latter.  Typical  conditions  of  interest  are  values  of  p  between  1-30  atm  and  values  of  T° 
between  1,000-1,800  K.  This  range  of  values  for  p  and  T°  was  used  to  guide  the  choice  of 
reactions  in  the  detailed  and  skeletal  mechanisms  and  the  initial  conditions  for  the  numerical 
problem  described. 

Computations  were  performed  using  the  SENKIN  computer  code  written  at  Sandia  National 
Laboratories  [12].  A  minor  modification  was  made  to  the  code  to  allow  input  of  the  functional 
form  for  pressure-dependent  rate  constants  (presented  earlier)  used  by  Tsang  and  Herron  [10].  A 
number  of  chemical  details  of  the  solutions  from  SENKIN  were  investigated  using  a 
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postprocessing  code  written  at  the  U.S,  Army  Research  Laboratory  (ARL)  [13].  Heat-release 
rate,  chemical  pathways,  and  logarithmically  normalized  sensitivities  of  the  calculated  species 
and  temperature  profiles  to  the  rate  constants  were  obtained  in  this  manner.  The  temperature 
sensitivities  were  quite  useful  because  these  qualitatively  reflect  the  sensitivity  of  the  calculated 
Tig  to  the  various  rate  coefficients  (information  not  directly  available).  The  pathway  and 
sensitivity  analyses  were  used  to  identify  those  reactions  that  have  a  significant  influence  on  the 
solution,  indicating  which  must  be  retained  in  the  skeletal  mechanisms. 

The  approximations  introduced  in  the  numerical  model  of  the  structure  of  the  dark  zone  are 
described  in  detail  elsewhere  [5, 9];  therefore,  only  a  brief  description  is  given  here.  The  system 
is  considered  to  be  spatially  homogeneous,  adiabatic,  and  isobaric.  Calculations  were  performed 
with  the  initial  mole  fractions  of  NO,  CO,  Hj,  Nj,  HjO,  and  COj  set  equal  to  0.24, 0.33, 0.08, 
0.04, 0.20,  and  0.10,  respectively,  for  double-base  propellents.  For  nitramine  propellents,  the 
initial  mole  fractions  of  NO,  CO,  H2,  Nj,  HjO,  COj,  NjO,  and  HCN  were  set  equal  to  0.13, 0.22, 
0.07, 0.06, 0.20, 0.09, 0.02,  and  0.22,  respectively.  The  results  presented  in  Table  2  demonstrate 
that  Tjg  calculated  using  the  skeletal  mechanism  agrees  very  well  with  the  values  obtained  from 
calculations  using  the  detailed  mechanism  [8]  for  the  double-base  propellents.  Values  of  Xjg 
calculated  for  the  nitramine  propellents  using  the  skeletal  mechanism  are  in  good  agreement  with 
those  of  the  detailed  mechanism  [8].  (There  are  differences  with  those  in  reference  8  due  to 
minor  revisions  of  the  detailed  mechanism). 

3.  REDUCED  CHEMICAL-KINETIC  MECHANISMS 

Procedures  for  deriving  reduced  chemical-kinetic  mechanisms  are  described  in  detail 
elsewhere  [14, 15].  The  skeletal  mechanisms  were  reduced  by  introducing  steady-state 
approximations  for  a  number  of  species.  For  any  species,  the  time  derivative  of  its  concentration 
is  the  difference  between  production  and  consumption  rates  for  that  species.  A  steady-state 
approximation  is  justified  if  either  the  species’  production  or  consumption  rate,  which  of  course 
are  nearly  equal  at  steady-state,  is  much  larger  than  the  time  derivative  of  its  concentration  [16]. 
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Table  2.  Comparison  of  the  Values  of  the  Ignition  Delay  Times  (Xjp  Calculated  Using 
the  Various  Mechanisms  for  Double-Base  Propellents 


Initial  Conditions 
(atm;  K) 

Detailed  Mechanism 
(41  sp.,  190  r.) 

(s) 

Skeletal 
Mechanism 
(15  sp.,  22  r.) 

(s) 

Reduced 
Mechanism 
(3  steps) 

(s) 

1;  1,000 

52.29 

51.47 

50.49 

1;  1,400 

1.035 

1.046 

1.042 

1;  1,800 

3.738  •  10-2 

3.843  •  10-2 

3.684  •  10-2 

5;  1,000 

8.297 

8.062 

7.769 

5;  1,400 

1.142  •  10-* 

1.147  •  10-* 

1.142  •  10-* 

5;  1,800 

7.517  •  10-2 

7.774  •  10-2 

7.624- 10-2 

10;  1,000 

3.994 

3.887 

3.725 

10;  1,400 

4.383  •  10-2 

4.375  •  10-2 

4.346-10-2 

10;  1,800 

3.448  •  10-2 

3.571  •  10-2 

3.503  - 10-2 

30;  1,000 

1.274 

1.255 

1.194 

30;  1,400 

9.90  •  10-2 

9.75  •  10-2 

9.63  - 10-2 

30;  1,800 

9.31  -lO-^ 

9.64- 10-^ 

9.49  - 10-^ 

Note:  The  initial  mole  fractions  of  NO,  CO,  H2,  N2,  H2O,  and  CO2  are  set  equal  to  0.24, 0.33, 0.08, 0.04, 0.20, 
and  0.10,  respectively. 

For  the  reduced  mechanisms,  the  concentrations  of  the  steady-state  species  were  calculated 
using  the  nonlinear  algebraic  equations  resulting  from  the  differential  rate  equations  for  these 
species  by  making  the  steady-state  assumptions.  That  is,  algebraic  equations  are  obtained  simply 
by  setting  these  species’  concentration  time  derivatives  to  zero.  The  steady-state  species’ 
concentrations  are  then  specified  in  terms  of  the  other  species’  concentrations.  However,  the 
initial  steady-state  species’  concentrations  are  no  longer  arbitrary.  For  example,  setting  them  to 
zero  is  internally  inconsistent.  This  causes  numerical  difficulties,  hi  the  calculations  with  the 
reduced  mechanisms,  the  steady-state  species’  initial  concentrations  were  therefore  first 
calculated  numerically  using  the  algebraic  equations.  For  all  cases  investigated,  the  initial 
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steady-state  species’  concentrations  were  found  to  be  negligibly  small  and  have  no  significant 
influence  on  the  calculated  Tjg. 

3.1  Double-Base  Propellents.  In  Figure  1,  the  ratios  of  the  time  derivatives  of  concentration 
to  production  rates  for  N,  NH,  HNO,  HONO,  H,  OH,  and  O  are  plotted  as  a  function  of  time  (t). 
The  calculations  were  performed  using  the  skeletal  mechanism  for  double-base  propellents  for 
p  =  5  atm  and  T®  =  1,400  K.  Figure  1  reveals  that  at  all  conditions,  except  those  close  to  t  =  0 
and  t  =  Tig,  the  value  of  this  ratio  is  very  small.  Therefore,  it  is  reasonable  to  introduce  steady- 
state  approximations  for  these  species.  In  Figure  2,  the  ratios  of  the  time  derivatives  of 
concentration  to  production  rates  for  NjO  and  NOj  are  plotted  as  a  function  of  time.  These 
results  were  also  obtained  using  the  skeletal  mechanism  for  double-base  propellents  for  p  =  5  atm 
and  T®  =  1 ,400  K.  Comparison  of  the  results  in  Figure  2  with  the  results  in  Figure  1  indicates 
that  the  concentration  of  NjO  is  not  in  steady-state,  but  it  is  still  reasonable  to  introduce  a  steady- 
state  approximation  for  NOj.  In  view  of  the  results  presented  in  Figures  1  and  2,  steady-state 
approximations  are  introduced  for  N,  NH,  HNO,  HONO,  H,  OH,  O,  and  NOj.  The  steady-state 
relations  provide  eight  nonlinear  algebraic  equations  used  to  calculate  the  concentration  of  the 
steady-state  species  in  terms  of  the  concentrations  of  the  remaining  species.  The  source  terms  of 
the  species  that  are  not  in  steady  state  are  rearranged  and  collected  in  three  groups,  such  that  each 
group  represents  the  rate  of  one  of  the  three  overall  steps  of  the  reduced  mechanism,  which  is 
written  as 


2NO  +  H2--N2O-1-H2O,  I 
N20-hH2-^N2  +  H20,  n 
and 

CO  +  H2O  ^  CO2 + H2.  in 

The  rates  of  the  overall  steps  w^,  k  =  I,  E,  III  expressed  in  terms  of  the  rates  of  the  elementary 
reactions  are 
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and 


=  Wj  + 

+  Wjj 

+ 

1 

Wn 

=  Wj  - 

W7 

W22, 

> 


w 


m 


=  w. 


(1) 


where  w„  is  the  net  reaction  rate  of  the  elementary  reaction  n  shown  in  Table  1. 


Figure  1.  The  ratio  of  the  difference  between  the  rates  of  production  and  the  rates  of 

consumption  to  the  rates  of  production  for  the  species  N.  NH.  HNO,  HONO,  H,  OH, 
and  O  plotted  as  a  function  of  time.  The  calculations  are  performed  using  the  skeletal 
mechanism  for  double-base  propellents  for  p  =  5  atm  and  T^ «=  1.400  K. 
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Figure  2.  The  ratio  of  the  difference  between  the  rates  of  production  and  the  rates  of 

consumption  to  the  rates  of  production  for  the  species  N^O  and  NO-,  plotted  as  a 
function  of  time.  The  calculations  are  performed  using  the  skeletal  mechanism  for 
double-base  propellents  for  p  =  5  atin  and  T”  =  1.400  K. 

3.2  Nitramine  Propellents.  Reduced  chemical-kinetic  mechanisms  for  nitramine  propellents 
were  obtained  using  the  same  procedure  employed  for  double-base  propellents.  The  rates  of 
production  and  consumption  of  various  species  were  calculated  using  the  skeletal  mechanism  for 
nitramine  propellents.  The  results  of  these  calculations  demonstrate  that  it  is  reasonable  to 
introduce  steady-state  approximations  for  NH,  HNO,  HNC,  NH2,  NOj,  OH,  and  O.  The  steady- 
state  assumptions  provide  seven  nonlinear  algebraic  equations  that  can  be  used  to  calculate  the 
concentrations  of  the  steady-state  species  in  terms  of  the  concentration  of  the  remaining  species. 
The  reduced  mechanism  has  six  overall  steps,  which  can  be  written  as 

2NO  +  2H2  -- N20  +  2H-hH20,  I 

H20-f-HCN  --  H2-HHNCO,  H 

N2O  +  2H2  N2  +  2H-I-H20,  m 
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CO  +  HjO  --  CO2  +  H2, 

2H  +  M  --  H2  +  M, 

NO  +  H  +  HNCO^  N2O  +  CO  +  H2. 


IV 

V 

VI 


The  rates  of  the  overall  steps  w,j,  k  =  I-VI  expressed  in  terms  of  the  rates  of  the  elementary 
reactions  are 


and 


— 

W2  Wg  +  Wji  , 

Wn 

W25  W33, 

Wm 

= 

Wj  +  W,  +  W22  +  W28  +  W32 , 

Wiv 

=: 

W3  +  W24  +  W32, 

Wy 

= 

W7  +  W12  +  W24  +  W28  W32, 

Wyi  =  W25 


+  W29  +  w. 


31 


(2) 


To  further  simplify  the  mechanism,  steady-state  approximations  were  introduced  for  the 
species  H  and  HNCO  yielding  four  overall  steps,  which  can  be  written  as 


2NO  +  H2  --  N20+H20, 

r 

NO  +  H2O  +  HCN  --  NjO  +  CO  +  1.5  H2, 

n' 

N2O  +  H2 --  N2  +  H2O, 

m' 

CO  +  H2O  --  CO2  +  H2. 

IV' 

The  rates  of  the  overall  steps  of  the  four-step  mechanism  w^,  k  =  F-IV'  are  identical  to  the  rates 
for  the  overall  steps  I-IV  in  the  six-step  mechanism  shown  in  equation  2, 
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4.  RESULTS  AND  DISCUSSION 


Calculations  with  the  reduced  mechanisms  were  performed  by  multiplying  time  derivative 
terms  in  the  differential  equations  for  the  steady-state  species  by  a  small  number.  Li  the  limiting 
case  where  the  value  of  this  number  is  zero,  the  differential  equations  are  transformed  to 
algebraic  equations.  The  value  of  the  multiplier  was  chosen  such  that  any  further  decrease  in  its 
value  had  no  influence  on  Xjg. 

In  Figure  3,  the  profiles  of  temperature  (T)  and  the  mole  fractions  of  NO,  Nj,  CO,  and  COj 
are  plotted  as  a  function  of  time.  The  calculations  were  performed  using  the  skeletal  mechanism 
and  the  reduced  three-step  mechanism  for  double-base  propellents  for  p  =  5  atm  and  T°  = 

1,400  K.  The  profiles  from  the  reduced  mechanism  agree  closely  with  those  calculated  using  the 
skeletal  mechanism,  except  for  times  close  to  ignition.  The  maximum  value  of  temperature  and 
the  final  values  of  the  mole  fractions  of  various  species  after  ignition  calculated  with  the  reduced 
mechanism  differ  from  those  calculated  using  the  skeletal  mechanism.  These  differences  occur 
because  the  concentrations  of  the  steady-state  species  are  not  included  in  the  atom-conservation 
equations  in  calculations  with  the  reduced  mechanism.  The  concentrations  of  the  steady-state 
species  represent  excess  mass  in  the  system.  Therefore,  in  the  calculations  with  the  reduced 
mechanism,  the  mass  fractions  of  the  various  atoms  change  with  time  and  the  final  mass  fractions 
of  these  atoms  are  different  from  their  initial  values.  However,  these  differences  become 
noticeable  only  at  times  close  to  Xjg.  In  fact,  calculations  prove  that  the  differences  between  the 
initial  mass  fractions  of  the  atoms  and  their  mass  fractions  at  times  approximately  90%  of  Xjg  are 
very  small. 

Figure  4  shows  profiles  of  temperature  (T)  and  the  mole  fractions  of  NO  and  N2  plotted  as  a 
function  of  time.  Figure  5  shows  profiles  of  T  and  mole  fractions  of  HjO  and  HCN.  These 
profiles  are  calculated  using  the  skeletal  mechanism  and  the  reduced  six-step  and  four-step 
mechanisms  for  nitramine  propellents  for  p  =  10  atm  and  T°  =  1,400  K.  The  profiles  in  Figures  4 
and  5  calculated  using  the  reduced  mechanisms  agree  closely  with  those  calculated  using  the 
skeletal  mechanism,  except  close  to  ignition. 
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Figures.  Profiles  for  temperature  fT)  and  mole  fractions  of  NO.  CO.  and  CO-,  plotted  as  a 

function  of  time.  The  calculations  are  performed  using  the  skeletal  mechanism  and  the 
reduced  three-step  mechanism  for  double-base  propellents  for  p  =  5  atm  and  T*’  = 


Table  2  shows  that  the  values  of  the  Xjg  calculated  using  the  reduced  three-step  mechanism  for 
double-base  propellents  agree  well  with  those  calculated  using  the  skeletal  mechanism.  The 
differences  in  the  values  of  Xjg  are  less  than  10%,  with  Tjg  from  the  reduced  mechanism  smaller 
than  that  obtained  using  the  skeletal  mechanism.  If  a  steady-state  approximation  is  introduced 
for  NjO,  a  reduced  two-step  mechanism  is  obtained.  The  value  of  Xjg  calculated  using  this 
two-step  mechanism  at  p  =  1  atm  and  T®  =  1,000  K  is  found  to  be  50%  lower  than  calculated 
with  the  three-step  mechanism.  This  result  illustrates  that  a  steady-state  approximation  for  NjO 
is  not  justified  (although  at  other  conditions  the  agreement  is  better). 


Table  3  compares  the  values  of  the  Xjg  calculated  using  the  six-step  and  four-step  mechanisms 
for  nitramine  propellents  with  those  calculated  using  the  skeletal  mechanism.  The  values  of  Xjg 


Mole  Fraction 


Figure  4.  Profiles  of  temperature  (T’)  and  mole  fractions  of  NO  and  N-,  plotted  as  a  function  of 
•  time.  The  calculations  are  performed  using  the  skeletal  mechanism  and  the  reduced 
six-step  and  four-step  mechanism  for  nitramine  propellents  for  p  =  10  atm  and  T”  = 
1.400  K. 


Figures.  Profiles  of  temperature  (T>  and  mole  fractions  of  H-,0  and  HCN  plotted  as  a  function 
of  time.  The  calculations  are  performed  using  the  skeletal  mechanism  and  the  reduced 
six-step  and  four-step  mechanism  for  nitramine  propellents  for  p  =  10  atm  and  T”  = 


1.400  K. 
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Table  3.  Comparison  of  the  Values  of  the  Ignition  Delay  Times  (tjg)- Calculated  Using  the 
Various  Mechanisms  for  Nitramine  Propellents 


hiitial 
Conditions 
(atm;  K) 

Detailed 
Mechanism 
(41  sp.,  190  r.) 

(s) 

Skeletal 
Mechanism 
(17  sp.,  23  r.) 

(s) 

Reduced 
Mechanism 
(6  steps) 

(s) 

Reduced 
Mechanism 
(4  steps) 

(s) 

1;1,000 

17.20 

17.34 

16.92 

16.60 

1;  1,400 

2.062  •  10-2 

1.995  •  10-2 

1.790  •  10-2 

1.322  •  10-2 

1;  1,800 

3.056  •  lO-'* 

2.602  •  lO"* 

2.160-  lO-'^ 

1.319-  lO-'^ 

10;  1,000 

2.754 

2.265 

2.185 

2.171 

10;  1,400 

5.723  •  10-^ 

6.162  •  10-2 

5.632  •  10-2 

4.388  - 10-2 

10;  1,800 

7.573  •  10-^ 

6.669  •  10-2 

5.590  •  10-2 

3.066  - 10-2 

30;  1,000 

1.241 

1.145 

1.099 

1.103 

30;  1,400 

3.122  •  10-2 

3.283  •  10-2 

3.037  •  10-2 

2.492  - 10-2 

30;  1,800 

4.903  •  10-2 

4.474  •  10-2 

3.803  •  10-2 

2.125-10:2  . 

Note:  The  initial  mole  fractions  of  NO,  CO,  N2,  H2O,  CO2,  N2O,  and  HCN  are  set  equal  to  0.13, 0.22, 0.07, 

0.06, 0.20, 0.09, 0.02,  and  0.22,  respectively. 

calculated  using  the  reduced  six-step  mechanism  agree  very  well  with  those  calculated  using  the 
skeletal  mechanism.  The  differences  are  less  than  18%.  However,  the  maximum  differences  in 
the  values  of  tjg  calculated  using  the  reduced  four-step  mechanism  and  the  skeletal  mechanism 
are  about  50%.  The  large  differences  between  the  values  of  Xjg  calculated  using  the  reduced 
four-step  mechanism  and  the  skeletal  mechanism  are  attributed  to  inaccuracies  introduced  by  the 
steady-state  approximations  for  H  and  HNCO. 

Pyrolysis  experiments  have  demonstrated  that  lower  temperatures  favor  production  of  CHjO 
and  NjO,  while  higher  temperatures  cause  more  HCN  and  NOj  to  be  formed  for  nitramines  and 
nitramine  propellents  [17, 18].  Clearly,  the  former  intermediates  reacting  near  the  propellent 
surface  are  expected  to  favor  formation  of  Nj  and  NjO  at  the  leading  edge  of  the  dark  zone,  while 
the  latter  favor  NO.  Because  a  higher  initial  dark-zone  temperature  would  be  expected  to 
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produce  a  steeper  temperature  gradient  near  the  propellent  surface,  temperatures  near  the  high 
end  of  the  range  studied  might  lead  to  initial  dark-zone  mixtures  that  contain  little  NjO. 
Therefore,  it  was  deemed  pmdent  to  ensure  that  if  future  experimental  data  exhibiting  this  trend 
become  available,  the  reduced  mechanisms  would  still  properly  predict  the  chemical  behavior  for 
such  mixtures.  For  this  reason,  Tjg  calculations  for  an  initial  temperature  of  1,800  K  were 
repeated  with  NjO  removed  from  the  initial  mixture.  The  results  are  shown  in  Table  4.  The 
differences  in  the  values  of  Tig  calculated  using  the  skeletal  mechanism  and  the  six-step  reduced 
mechanism  are  less  than  17%,  and  the  maximum  difference  in  the  values  of  Xjg  calculated  using 
the  skeletal  mechanism  and  the  four-step  reduced  mechanism  is  about  60%. 


Table  4.  Comparison  of  the  Values  of  the  Ignition  Delay  Times  (Xig)  Calculated  Using  the 
Various  Mechanisms  for  Nitramine  Propellents 


Initial 
Conditions 
(atm;  K) 

Detailed 
Mechanism 
(41  sp.,  190  r.) 

(s) 

Skeletal 
Mechanism 
(17  sp.,  23  r.) 

(s) 

Reduced 
Mechanism 
(6  steps) 

(s) 

Reduced 
Mechanism 
(4  steps) 

(s) 

1;  1,800 

2.035  •  10-^ 

1.446  •  10'^ 

1.204  •  10-^ 

5.738  •  10-^ 

10;  1,800 

5.858  •  10'^ 

5.233  •  10-^ 

4.452  •  lO"* 

3.026  •  10-^ 

30;  1,800 

3.270  •  10'^ 

3.090  •  10-^ 

2.671  •  lO"* 

1.944*  10-^ 

Note:  The  initial  mole  fractions  of  NO,  CO,  H2,  N2,  H2O,  CO2,  and  HCN  are  set  equal  to  0.133, 0.224, 0.071, 
0.061, 0.204, 0.091,  and  0.224,  respectively. 


The  data  in  Tables  3  and  4  highlight  the  importance  of  the  trace  of  NjO  in  the  initial  mixture. 
Comparison  of  results  in  these  tables  shows  that,  in  general,  the  presence  of  the  N2O  decreases  Xjg 
by  nearly  a  factor  of  2  at  1,800  K.  The  effect  is  even  greater  at  lower  temperatures  (not  shown). 
Addition  of  NjO  to  the  mixtures  results  in  an  increase  in  the  rate  of  the  reaction  NjO  +  M  =  Nj 
+  O  +  M  (Table  1,  elementary  step  1),  an  important  radical  source. 
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5.  SUMMARY  AND  CONCLUSIONS 


Skeletal  and  reduced  chemical-kinetic  mechanisms  are  derived  for  predicting  the  structure  of 
the  dark  zones  formed  over  the  burning  surface  of  typical  double-base  and  nitramine  propellents. 
In  the  calculations  with  the  reduced  mechanisms,  the  final  mass  fractions  of  the  atoms  after 
ignition  differ  from  their  initial  values  because  the  concentrations  of  the  steady-state  species  are 
not  included  in  the  atom  conservation  equations.  Therefore,  the  reduced  mechanisms  are  suitable 
for  calculating  the  values  of  tjg,  but  in  some  applications,  the  errors  introduced  in  the  final 
temperature  and  species’  concentrations  may  not  be  acceptable.  Attempts  to  improve  the 
accuracy  of  the  predictions  of  the  final  concentrations  of  these  species  in  the  reduced 
mechanisms  are  in  progress.  It  is  hoped  that  these  skeletal  and  reduced  chemical-kinetic 
mechanisms  will  be  suitable  for  use  in  interior  ballistics  calculations.* 


*  Recently,  we  have  perfonned  further  updates  to  the  detailed  chemical  mechanism  from  which  the  skeletal  and 
reduced  mechanisms  are  derived.  We  have  discovered  that  Qiree  important  reactions  involving  two  additional 
species  were  overlooked  for  the  detailed  mechanism.  These  changes  have  pronounced  effects  at  low 
temperatures,  especially  at  high  pressures,  resulting  in  shortening  of  predicted  ignition  delays  by  up  to  a  factor  of 
3.  The  effects  are,  unfortunately,  negligible  in  regions  where  experimental  data  are  available  for  testing.  Reduced 
mechanisms  have  been  developed  from  the  newest  detailed  mechanism  in  a  way  entirely  analogous  to  that  used 
herein.  The  results  are  qualitatively  similar  in  terms  of  accuracy  compared  to  the  detailed  mechanism.  The 
number  of  reactions  in  both  skeletal  mechanisms  increases  by  three  and  the  steady-state  species  increase  by  two. 
These  new  results,  as  well  as  comparisons  to  our  previous  results,  wiU  be  presented  in  future  reports. 
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